Susceptibility of the QCD vacuum to CP-odd electromagnetic background fields 
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We investigate two flavor QCD in presence of CP-odd electromagnetic background fields and 
determine, by means of lattice QCD simulations, the induced effective 6 term to the first order in 
E ■ B. We employ a rooted staggered discretization and study lattice spacings down to 0.1 fm, 
with pion masses around 400 MeV. In order to deal with a positive measure, we consider purely 
imaginary electric fields and real magnetic fields, then exploiting analytic continuation. Our results 
are relevant to a description of the effective pseudoscalar QED-QCD interactions. 

PACS numbers: 12. 38. Aw, 11.15.Ha,12.38.Gc 



a; 
43- 



> 

(N 
(N 

p; 

ov 
o: 

(N 



X 



I. INTRODUCTION 

Quantum ChromoDynamics (QCD) may contain in- 
teractions which violate CP, the symmetry under charge 
conjugation and parity, corresponding to a term —iOQ 
in the Euclidean action, where 



Q = I d 4 x q(x) = 



d 4 x- 



64tt 



:G^(x)G« v (x) (1) 



is the topological charge operator, is the non-Abelian 
gauge field strength and G a = e^upoG®. However, ex- 
perimental upper bounds on the parameter 9 are quite 
stringent, \0\ < 10" 10 [H, @j. Nevertheless, 6 related ef- 
fects play an important role in strong interaction phe- 
nomenology and are generally linked to fluctuations of 
Q, which affect, through the axial anomaly, the balance 
of chirality. 

A significant interest has been raised recently by the 
possibility that local effective variations of 6, correspond- 
ing to topological charge fluctuations, may induce de- 
tectable phenomena in presence of magnetic fields as 
strong as those produced in the early phases of non- 
central heavy ion collisions, which can reach up to 10 15 
Tesla at LHC. According to the so-called chiral magnetic 
effect [3- 5] , the net unbalance of chirality induced by the 
topological background would lead, in presence of a mag- 
netic field strong enough to align the magnetic moments 
of quarks, to a net separation of electric charge along the 
field direction. 
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The physics of strong interactions in presence of elec- 
tromagnetic (e.m.) backgrounds has attracted much in- 
terest also in relation to the possible effects of the back- 
ground field on the QCD vacuum and on the QCD phase 
diagram, stimulating a number of model computations 
and, more recently, of lattice simulations. Such effects 
may be relevant in various contexts: magnetic fields of 
the order of 10 16 Tesla may have been produced at the 
cosmological electroweak phase transition [6] ; large mag- 
netic fields are also expected in compact astrophysical 
objects such as magnetars 

One aspect, emerging from recent lattice simulations 
with dynamical fermions [§l-[l3T] and from model stud- 
ies [14l4l6l |. is that e.m. background fields, even if directly 
coupled only to charged particles, may have a significant 
influence, via quark loop effects, also on the gluonic sec- 
tor. In this study, in an attempt to better clarify such is- 
sue, we will investigate how the explicit breaking of some 
symmetry by the e.m. background field propagates to 
gluon fields, considering the particular case of CP sym- 
metry. 

To be more specific, let us consider QCD in pres- 
ence of a constant and uniform e.m. field such that 
F^vF^ oc E ■ B ^ 0: such background is expected to 

induce an effective CP-violating interaction in the gluon 

2 

sector, ^effgf^a G^ v (x)G't v {x) . 8 c s must be an odd func- 
tion of E ■ B and, at the lowest order, we can write 



eff ~ XcP e 2 E-B + 0{{E-B) 



(2) 
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where xcp is a sort of susceptibility of the QCD vacuum 
to CP-breaking e.m. fields. 

The effect that we want to study is in some sense com- 
plementary to the chiral magnetic effect, where a CP- 
violating non-Abelian background gives rise to charge 
separation, hence to an electric field, parallel to a back- 
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ground magnetic field: in that case CP violation propa- 
gates from the gluonic to the e.m. sector, i.e. opposite 
to what we investigate here. In fact, xcp is directly re- 
lated to the strength of the effective pseudoscalar QED- 
QCD interaction, xcp q(x) e 2 E ■ B [r7 Hl9|; i n particu- 
lar, to connect with the notation of Ref. we have 
Xcp — K /2, with k defined as in Eq. (5) of Ref. (l9| . 

The purpose of this study is to furnish a first deter- 
mination of xcp based on lattice QCD simulations. To 
that aim we will do simulations of QCD in presence of 
uniform e.m. background fields such that E ■ B ^ 0, and 
then we will determine the corresponding induced 9 c s by 
studying the topological charge distribution. 



II. THE METHOD 

An e.m. background field enters the QCD lagrangean 
by modifying the covariant derivative of quarks, = 
+ igA®T a + iqA^, where is the e.m. gauge po- 
tential and q is the quark electric charge. In the simplest 
discretization on a cubic lattice, that is implemented by 
adding proper U (1) phases M M (n) to the usual SU(3) par- 
allel transports, U^n) — > w M (n)[/ M (n), where n is a lat- 
tice site. 

A constant and uniform e.m. field with a single non 
vanishing component = F can be realized by a po- 
tential A v = FXfi and A p = for p ^ v. In presence of 
periodic boundary conditions (b.c), which is the usual 
choice in lattice simulations to minimize finite size effects, 
F must be integer multiple of a minimum quantum 



/ = 2ir/(qa 2 L fi L l/ ) , 



(3) 



and proper b.c. must be chosen for fermions, to preserve 
gauge invariance [20| . The corresponding U(l) lattice 
links are 



ul q) (n) 



,('/) 



-ia qL^Frii, 



(4) 



and u p (n) = 1 otherwise, where L M is the number of 
lattice sites in the /i direction. In presence of various 
non- vanishing components of F^, the definition above 
generalizes by simply writing the phase of each U(l) link 
as the sum of the contributions from the various non- 
vanishing components of F^ LU . 

In the following we shall consider two flavor QCD with 
fermions discretized in the standard rooted staggered for- 
mulation. In the corresponding functional integral, each 
quark is described by the fourth root of the fermion ma- 
trix determinant: 



Z = J VUe- s ° detM*[U,q u ]detMz[U,q d } 



(5) 



1,3 



^ 4 

- u<«\i-v) U}(i-i>)5 id+i> ) . (6) 



T>U is the functional integration over the non-Abelian 
gauge link variables, Sq is the plaquette action, i and j 
refer to lattice sites and rj u (i) are the staggered phases. 
The choice for the quark charges is standard, q u — 2|e|/3 
and qd = — |e|/3, leading to a quantization in units of 
/ = 6ir/ (|e|a 2 L^L^) for each component of F^ v . 

The addition of pure U(l) phases to the standard 
SU(3) link variables leaves the spectrum of M — mid 
purely imaginary and symmetric under conjugation: that 
guarantees detM > 0, hence the feasibility of numerical 
simulations. However, it is easy to realize that -Fbz ^ 0, as 
defined above, corresponds, when continued to Minkovski 
space, to a purely imaginary electric field pll [22l] . 

In order to have a real electric field, one should intro- 
duce imaginary components for the gauge potential in 
Euclidean space, but that would take the m m variables 
out of the U{\) group and make the quark determinant 
complex, thus introducing a sign problem which hinders 
numerical simulations. On the other hand, this is ex- 
pected in view of the phenomenon that we want to ex- 
plore since, as it is well known, also the formulation of 
QCD at real non-vanishing 9 suffers from a sign problem. 

The usual attitude, taken e.g. in lattice studies of 
the electric polarizabilities of hadrons, is to keep the 
electric background field purely imaginary, to avoid the 
sign problem, then exploiting analytic continuation. We 
shall follow the same approach, making use, in Euclidean 
space, of the same definition for all components of F pv : 
that will correspond to real magnetic fields B and imag- 
inary electric fields E = i Ej. As a consequence, we 
expect to produce a purely imaginary effective parame- 
ter e ff = i9i e g. On the other hand, working with an 
imaginary 9 is also one of the possible approaches to the 
study of ^-dependence in QCD (2^3. 

The presence of an imaginary 9\ adds a factor 
exp(9jQ) to the probability distribution of gauge fields. 
That will shift the distribution of the topological charge 
by an amount which, at the linear order in 9i, is governed 
by the topological susceptibility x at 9j = 0: 



(Q)e I ^V X di = (Q 2 )e=oO I 



(7) 



where V is the spacetime volume. That gives us the 
opportunity of determining the effective 0/ e ff produced 
by a given e.m. field as 



7/cff 



(Q)(Ei,B) 

{Q 2 )o 



0((E! ■ B) 



(8) 



where by (-)o we mean the average taken at zero e.m. 
field. Corrections to Eq. will be negligible at least 
in the region of small 9i c g which is relevant to Eq. ©. 
In the following we will show that the so defined 9i e g 
is indeed an odd function of Ej ■ B alone and, assuming 
that the theory is analytic at least for small enough back- 
ground fields, we shall determine the susceptibility xcp 
defined in Eq. ([2]) from the small field behavior of 9i e s- 
For the determination of Q on gauge configurations, 
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5000 10000 
RHMC trajectories 

FIG. 1: Monte-Carlo history of Q for ~ 480 MeV and a ~ 
0.12 fm, on a 16 4 lattice, for two different e.m. backgrounds. 



we shall adopt a standard discretized gluonic definition, 
-1 ±4 



fivpa— ±1 



measured after cooling [2g, 129J, i.e. minimization of 
the pure gauge action to eliminate ultraviolet (UV) ar- 



: Ejui/per for positive directions and 
Such method is based on the as- 



tifacts; here e M „ p(T 

(pvpa — £( — fj,)upcr- 

sumption that a given configuration remains in the start- 
ing topological sector as the action is minimized; it is 
known to provide correct results, comparable with im- 
proved fermionic definition, as the continuum limit is 
approached, while maintaining a limited computational 
cost. 



III. RESULTS 

We have studied N f =2 QCD at T = (see Eq. ©), 
for a fixed pseudo-Goldstone pion mass ~ 480 MeV. 
Different lattice spacings have been explored by tuning 
the inverse gauge coupling j3 and am according to what 
reported in Ref. (30|: a ~ 0.30 fm (/3 = 5.35, am = 
0.075), a ~ 0.20 fm {(3 = 5.453, am = 0.02627), 
a ~ 0.15 fm (0 = 5.527, am = 0.0146), a ~ 0.12 
fm (j8 = 5.5829, am = 0.01048) and a ~ 0.10 fm 
(/9 = 5.6286, am = 0.008318). Different lattice volumes 
have been explored to check for finite size corrections (see 
Pig.®. 

Numerical simulations have been performed by a stan- 
dard Rational Hybrid Monte-Carlo (RHMC) algorithm 
on GPU cards, making use of the code illustrated in 
Ref. [HJ. Typical statistics have been of the order of 
10K molecular dynamics time units for each run. The 
number of (straight) cooling sweeps n coo i has varied from 
30 for the smallest to 60 for the largest spacing adopted, 



FIG. 2: Distribution of the topological charge after 30 cooling 
steps. Data are the same as for Fig. [T] The factor a (see text) 
is 0.947 in this case. 



and it has been chosen in a window where results show 
no significant dependence on n coo i. After cooling, the 
charge is first divided by a constant factor a close to 
one, so that its distribution gets peaked around integer 
values (see e.g. Fig. [2]), then it is approximated to the 
closest integer to finally get Q. Autocorrelation times for 
Q, which are known to grow critically when approaching 
the continuum limit, have reached up to O(100) RHMC 
trajectories for a ~ 0.1 fm. 

As a first illustration of our results, in Fig. [T]we show 
the Monte-Carlo history of Q for two numerical simula- 
tions performed respectively at zero and non-zero e.m. 
background field, on a 16 4 lattice with a ~ 0.12 fm and 
ttItt ~ 480 MeV. The non-zero e.m. field corresponds to 
Ei = B = 3/z, where / is defined in Eq. ([3]), corre- 
sponding to e 2 Ei ■ B ~ 0.37 GeV 4 . We see that, while in 
absence of the e.m. background Q fluctuates around zero, 
as expected by CP-invariance, fluctuations are shifted 
towards positive values as Ei ■ B ^ 0. This is more 
clearly visible in Fig. [2] where we plot the correspond- 
ing distributions of Q; in this particular case we obtain 
{Q)(E I ,B)/{Q 2 ) = 1.46(14). 

In order to better investigate the dependence of 
(Q)(Ej , B) on the background field values, in Fig. [3] we 
show (Q)(E I ,B)/(Q 2 ) for a 16 4 lattice with a ~ 0.3 fm 
and m^ ~ 480 MeV. Data are obtained for a variety of 
combinations of Ej and B, mostly taken parallel to the 
z axis, and then plotted versus Ei ■ B. The fact that all 
data fall on the same curve, even when Ej and B are not 
parallel, is a nice demonstration that, within errors, #j e ff 
is indeed a function of Ej ■ B alone, as expected; we have 
verified that explicitly by taking different combinations 
of the fields having exactly the same or opposite values 
for Ei ■ B, thus checking also that 8i c s is odd in Ej ■ B. 
The dependence is linear in Ej ■ B for small fields, then 
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FIG. 3: (Q)(Ei,B)/{Q 2 ) for various (Ei, B) on a 16 4 lat- 
tice for rriir ~ 480 MeV and a ~ 0.3 fm. Open circles corre- 
sponds to the following combinations of Ei-B/f 2 : (2, —8, —2)- 
(8, 1, 4), (1, 2, 3) ■ (3, 2, 1), (1, 3, 4) ■ (2, 2, 3), (5, 5, 5) ■ (1, 2, 3), 
(6, 6, 6) ■ (6, 6, 6), and (8, 8, 8) • (8, 8, 8). The dashed line cor- 
responds to a best fit to Eq. (110[) . 
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FIG. 4: B)/(Q 2 ) as a function of Ej ■ B for ~ 

480 MeV and different spacings o and lattice volumes L 4 . 
The continuous line corresponds to a best fit to Eq. (| 10p at 
the smallest value of a, the dashed line is the corresponding 
slope at Ei ■ B — 0. 



showing signs of saturation for larger fields, as common 
to many systems with a linear response to external stim- 
ulation. We have found that all data can be nicely fitted 
by a function 

(Q)(E Il B)/(Q 2 ) = a atan( ai S • B) ; (10) 

the best fit curve, corresponding to % 2 /d.o.f. = 0.74, is 
shown in Fig. [3] 

In order to discuss finite size and finite UV cutoff ef- 
fects, in Fig. H we show {Q)(Ei, B)/(Q 2 ) for m n ~ 480 
MeV and different spacings a and lattice volumes L A . 
(Q) and (Q 2 )o are both derivatives of the free energy 
with respect to 9, hence they are proportional to V and, 
apart from possible systematic effects, their ratio should 
be independent of V. From Fig.|4]we infer that finite size 
effects are not significant, even on the smallest volumes 
explored, corresponding to am n L ~ 4. 

The dependence on the UV cutoff instead seems sig- 
nificant until a < 0.15 fm. Apart from standard lattice 
artifacts related to the path integral discretization, ad- 
ditional systematic effects may be related to the method 
used to determine Q: if a is so coarse that part of the 
topological background, created by the influence of the 
e.m. field, lives close to the UV scale, then cooling is 
expected to destroy part of such background. However 
data obtained for a < 0.15 fm are in fairly good agree- 
ment with each other, within errors, especially in the 
small field region, which is the one relevant for the deter- 
mination of xcp ■ 

In order to determine xcp j we have performed best fits 
of the data in Fig.|4]to the function in Eq. (flOl) . in a range 
e 2 Ej-B < 0.6 GeV 4 , then considering its slope at Ej-B = 
and exploiting Eqs. ((2J and ((HJ). We have verified that, 
for each set, the slope is consistent, within errors, with 



a direct linear fit restricted to a narrow enough region of 
small Ei -B. We obtain, for the smallest lattice spacings, 
Xcp = 7.1(6) GeV" 4 for a = 0.15 fm, X cp = 6.9(5) 
GeV" 4 for a = 0.12 fm and xcp = 7.2(7) GeV~ 4 for 
a = 0.10 fm, so that we can quote xcp — (7± 1) GeV -4 
as a reasonable estimate for this value of the pion mass. 
Preliminary results obtained on a 16 4 lattice and for a ~ 
0.15 fm indicate instead xcp — 10(1) GeV -4 if ~ 
280 MeV, suggesting that xcp tends to increase when 
approaching the chiral limit. 



IV. DISCUSSION 

It is useful to compare our results with the phe - 
nomenological estimate of K = 2 xcp given in Rcf. 191, 
which is based on the effective couplings of the 77 and 
T]' mesons to two photons and to two gluons: xcp ~ 
0.73/(Tr 2 f 2 m 2 ,) ~ 3 GeV -4 : even considering the differ- 
ent systematics (the phenomenological estimate is based 
on a theory with 2+1 light flavors), our data suggest 
that the lattice QCD result for the effective pseudoscalar 
QED-QCD interaction is larger, even if of the same order 
of magnitude. 

Regarding the validity of analytic continuation from 
imaginary to real electric fields, we notice that, while a 
smooth behavior is expected as the imaginary electric 
field approaches zero, at least for non-zero quark masses, 
a non-zero and constant real electric field, however small, 
is instead expected to induce vacuum instabilities. On 
the other hand this is not true in presence of an infrared 
cutoff, i.e. if electric fields are limited in space. Therefore 
our result should be applicable to determine the local 
effective 9 parameter produced by smooth enough but 
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limited in space CP-violating e.m. fields. It would be 
interesting in the future to consider the case of smoothly 
varying fields explicitly. 

In order to give an idea of the magnitude of the effect 
that we have studied, our estimate predicts that two par- 
allel magnetic and electric fields, with eB ~ 1 GeV 2 and 
eE ~ 1 MeV 2 , would induce an effective 9 of the order 
of 1CT 5 . 

Finally, apart from repeating our determination with 
more physical quark masses and closer to the continuum 
limit, it would be interesting to extend our investiga- 
tion also to finite temperature, especially across and right 
above the deconfinement transition, where the determi- 
nation of the effective pseudoscalar QED-QCD interac- 



tion would be relevant also to the phenomenology of the 
chiral magnetic effect in heavy ion collisions. 
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